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We report on our ongoing project to calculate the leading hadronic contribution to the anomalous 
magnetic moment of the muon a^ LO using two dynamical flavours of non-perturbatively &{d) 
improved Wilson fermions. In this study, we changed the vacuum polarisation tensor to a 
combination of local and point-split currents which significantly reduces the numerical effort. 
Partially twisted boundary conditions allow us to improve the momentum resolution of the 
vacuum polarisation tensor and therefore the determination of the leading hadronic contribution 
to (g — 2)^. We also extended the range of ensembles to include a pion mass below 200 MeV 
which allows us to check the non-trivial chiral behaviour of a^ L °. 
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1. Introduction 

The anomalous magnetic moment of the muon = (g^ — 2)/2 currently shows a persistent 
discrepancy of arround 3 sigma between experiment and theoretical predictions which may be a 
hint to physics beyond the Standard Model [1, 2]. At the moment can be used to put exclusion 
limits on several BSM models, among many we mention the search for dark photons, on which the 
experimental groups at the KPH Institute in Mainz are involved [3]. The theoretical uncertainties 
of anomalous magnetic moment of the muon are dominated by hadronic contributions, from which 
the lowest order hadronic vacuum polarisation fl^ LO accounts for the largest uncertainty. At the 
moment the Standard Model prediction of is obtained by using the optical theorem to relate 
the e + e~ — > hadrons cross-section data to the vacuum polarisation or by using the measured spectral 
functions from hadronic T decays. We report on our ongoing project to determine the leading 
hadronic contribution to the anomalous magnetic moment of the muon from first principles using 
Lattice QCD. 

2. Lattice Setup 
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Table 1: Summary of simulation parameters. The pion masses and the scale have been taken 
from [13, 14, 15] and are partly still preliminary. 

We use non-perturbatively &{a) improved Wilson fermions with two dynamical flavours [4] 
and include a strange quark in a partially quenched framework. Our simulations were performed on 
a subset of the gauge configurations generated within the CLS collaboration [5]. The corresponding 
simulation parameters are listed in table 1, where the ensembles highlighted in red have been 
recently generated and added to the set analyzed in our earlier publication [6]. Similar studies have 
been performed in the quenched approximation [7, 8] and in the theory with two [9], three [10,11] 
and four dynamical flavours [12] . 

3. Vacuum Polarisation 

On the lattice in Euclidean space-time, the vacuum polarisation tensor is defined by a Fourier 
transformation of a current-current correlator 

U, v (q) = Zya^e^" 1 ' 2 ) (^(x)4(0)) , (3.1) 
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in which we use a combination of the local current J\{x) = q{x)y ll q{x) and conserved point- 
split current J^(x) as done in [11]. The conserved point-split current J^(x) can be derived using 
Noether's theorem and reads for Wilson fermions as 

W = \ (q(x+apL)(l + Yv)U+(x)q(x) - q(x)(l - Yfl )U^x)q(x + afl)) . (3.2) 

The Wick contraction of eq. (3.1) leads to connected and disconnected contributions. In this work 
we currently only consider the connected contributions, nevertheless two-flavour #PT can be used 
to estimate the remaining disconnected piece to be approximately —10% of the connected one [16]. 
The vacuum polarisation H(q 2 ) can be extracted using the following relation: 

n Mv(<?) = (<W 2 - q^qv) n(g 2 ). (3.3) 

The use of the local current in the definition of the vacuum polarisation tensor requires a renormal- 
ization of the operator. The relevant renormalization constant Zy was determined independently 
in [17]. Using the combination of local and conserved currents is computationally cheaper. Fig- 
ure 1 shows a comparison of our previous work [6] using two conserved point-split currents and 
the combination of local and conserved current on the same gauge configurations. The unphysical 
constant 11(0) turns out to be different in both determinations. Nevertheless the subtracted vac- 
uum polarisation which enters eq. (3.4) agrees very well within the given statistical precision. The 



TT(q 2 
1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




PRELIMINARY 



local - conserved 
conserved - conserved 



0* 



0.2 



0.4 0.6 
q 2 [GeV 2 ] 



0.8 



Figure 1: Comparison of the subtracted vacuum polarisation using two conserved point-split currents (red) 
and a combination of local and conserved point-split current (blue). The subtracted vacuum polarisation is 
shown for the F6 ensemble (J3 = 5.3, L = 3.0fm, m % = 324 Me V). 



lowest order hadronic contribution to the anomalous magnetic moment of the muon a^ h ° can be 
obtained by integrating 

dq 2 f(q 2 )tl(q 2 ), (3.4) 
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n) jo 

which is the product of the subtracted vacuum polarisation tl(q 2 ) = 471 (Tl(q 2 ) — 11(0)) and a func 
tion originating from QED [18, 19] 
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Because of the shape of f(q 2 ), the low momentum region of the vacuum polarisation gives the most 
significant contribution to a^ L °. To improve the momentum sampling we implemented twisted 
boundary conditions to the valence quarks [20] 

q(x + Lk) =e Wk q(x). (3.6) 

This modifies the momentum to ^H— £, where can be tuned to any real value. Using isospin 
symmetry, the connected contribution of the vacuum polarisation can be understood as flavour 
non-diagonal and partially twisted boundary conditions can be applied to evaluate it (see [16]). 

4. Determination of aS L0 



n(q 2 ) 




Figure 2: Left: The vacuum polarisation Tl(q 2 ) in the low-momentum region, computed on the G8 ensemble 
(/3 = 5.3, L = 4.0fm, m n = 190 Me V) using twisted and periodic boundary conditions. The green curve 
shows the result for a correlated fit to a Pade eq. (4.4). Right: The contribution to aJJ LO from the integral in 
eq. (3.4) separated by the different momentum regions. 

The leading order hadronic contribution to can be obtained by integrating eq. (3.4) numer- 
ically. For this a continuous description of the vacuum polarisation H(q 2 ) is required. We perform 
correlated least-square fits to the simulation data. Through varying the fit ansatz we are able to 
check for systematic uncertainties, for which we choose 

a) a vector dominance model including a single vector 

b) a vector dominance model with two vectors and one mass fixed to my as proposed in [1 1] 
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c) an ansatz motivated by a model for the hadronic cross section as used in [8] 

U(q 2 ) = a + b-log(q 2 + c 2 ) + d (4.3) 

Or + e ) 

d) a model independent Pade with various degrees as proposed in [6, 21, 22], e.g. 

For the fit ansatz we imply that perturbation theory is matched smoothly at large momentum q 2 > 
2GeV 2 , which allows us to reduce the number of fit parameters by one. Further details for the 
fitting procedure and the matching to perturbation theory can be found in [6]. All fits, except the 
vector dominance model fit eq. (4.1), agree within the statistical uncertainties. Figure 2 shows both 
the fit results for the vacuum polarisation for small momentum and the clear improvement in the 
momentum resolution achieved by partially twisted boundary conditions. The right panel shows 
the individual contributions to <3^ LO from different momentum ranges. It is important to note that 
without twisted boundary conditions only a single data point lies in the momentum range covered 
by regions I and II. 



5. Results 
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Figure 3: Left: The simulation results for a^ LO computed using two flavours, shown as function of nc\. 
The chiral extrapolation (blue curve) is performed using an ansatz motivated by chiral perturbation theory. 
Right: Equivalent results for a^° including a partially quenched strange quark. 

Repeating the fitting procedure for each of the ensembles listed in table 1 we obtain results for 
the chiral behavior a^ LO shown in figure 3. The data show a clear non-linear behaviour and a rise 
towards the physical pion mass. A chiral extrapolation is necessary to extract a^ LO at the physical 
point. We use an uncorrelated fit to an ansatz motivated by chiral perturbation theory, i.e. 

4 L °(»4) =A+S-m 2 +C-m 2 log(m 2 ). (5.1) 
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In addition we perform a linear fit to the most chiral points 

a* L0 (ml)=A + B-m 2 K (5.2) 

to check for systematic uncertainties in the chiral extrapolation. Figure 3 shows the corresponding 
fits as well as the uncertainties. We find that the #PT motivated fit describes the entire range of 
data points quite well. With the current level of accuracy, the three flavour extrapolation misses 
the phenomenological value for a^ h0 , which might be caused by statistical downward fluctuations 
of the most chiral points. To verify the chiral behaviour we will increase our statistics on these 
ensembles. 

6. Conclusions and Outlook 

The determination of the leading hadronic contribution to the anomalous magnetic moment of 
the muon using lattice QCD has received considerable attention in the last years, but still requires 
additional improvements to have an impact on phenomenology. The statistical uncertainty of the 
individual data points is currently of order 2 — 7%, where the largest uncertainties is obtained for 
the most chiral data points. Systematic uncertainties arising from chiral extrapolations and lat- 
tice artefacts will come on top for the final lattice estimate. Partially twisted boundary conditions 
have proven to improve the sampling of the momentum for the vacuum polarisation and thereby 
reduce the statistical and systematic uncertainties. The combination of conserved and local vector 
current reduces the numerical effort, which allows one to improve the accuracy of the determina- 
tion. Further details of our study will be published soon [23]. Once all systematic uncertainties 
are under control, we will include a dynamical strange and charm quark into our simulations. By 
combining the strategy described here and "all mode averaging" [24, 25] plus the approach in [26] 
to directly estimate 11(0), the connected part of the leading order hadronic contribution to (g — 2)^ 
can be computed at level of less than few percent. At that point the dominating uncertainty will 
come from the missing disconnected contribution. The application of the techniques in [27] looks 
promising in this respect. 
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